Non-equilibrium steady-states for interacting open systems: exact results 
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Under certain conditions we prove the existence of a steady-state transport regime for interacting 
mesoscopic systems coupled to reservoirs (leads). The partitioning and partition-free scenarios are 
treated on an equal footing. Our time-dependent scattering approach is exact and proves, among 
other things the independence of the steady-state quantities from the initial state of the sample. 
Closed formulas for the steady-state current amenable for perturbative calculations w.r.t. the inter- 
action strength are also derived. In the partitioning case we calculate the first order correction and 
recover the mean-field (Hartree-Fock) results. 
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INTRODUCTION 



The theoretical modeling of time-dependent transport 
has been an active area of research in the last few 
years^^^. Transient currents are calculated using the 
Keldysh formalism^, and the electron-electron interac- 
tion (EEI) effects are accounted for via time-dependent 
density functional theory (TDDFT) or many-body per- 
turbative (MBP) methods. But no matter which metod 
one uses and how simple the system is, one recovers an old 
and often avoided question in non-equilibrium transport: 
does there exist a non-equilibrium steady-state (NESS), 
and if yes, is it unique? 

Recently, Kurth et al.^ presented TDDFT simulations 
for a single-level quantum dot (QD) in the Coulomb 
blockade regime. By selecting a suitable exchange- 
correlation potential they noticed that the system does 
not evolve to a steady state but rather follows charg- 
ing/discharging cycles. Moreover, Myohanen et al.^ and 
Puig et al.^ emphasized that different approximation 
schemes for the interaction self-energy lead to different 
steady-states, i.e, in the long-time limit the numerical 
simulations lead to different values of the current. These 
recent findings cleary show that the crossover to a steady- 
state (if any) is a non-trivial aspect which is revealed only 
by a fully time-dependent formalism for open and inter- 
acting systems. We remind here that following Ref.^ the 
Keldysh formalism was extensively used to compute the 
stationary currents by assuming i) that such a steady 
state is achieved and ii) that the interaction strength is 
rather small such that a perturbative approach makes 
sense. The first assumption implies that it is sufficient to 
work directly with the Fourier transforms of the Green 
functions. The second assumption allows one to exploit 
diagrammatic techniques and conserving approximations 
for the interaction self-energy^. Given the results men- 
tioned above one faces three questions: 1) How legitimate 
is it to take for granted the steady-state quantities of the 
Keldysh formalism especially in the presence of electron- 
electron interaction ? 2) Is the result of Kurth et al.^ 



universal? 3) Is it possible to establish the existence of a 
stationary current in the long-time limit before selecting a 
given approximation scheme for the explicit calculations 
of the interaction effects? 

In this note we prove that an interacting sample evolves 
to a NESS provided: i) all single-particle eigenstates 
of the isolated sample become resonances (with positive 
width) when the leads are coupled and ii) the interaction 
strength is small enough to ensure convergence of a cer- 
tain perturbation expansion. These conditions are met 
for large quantum dots coupled to broad leads. Condition 
i) is also fulfilled if the bias applied on the leads covers 
the entire spectrum of the sample (wide-band limit). We 
follow the scattering approach to the NESS of open quan- 
tum systems advocated by Ruelle^^ and implemented, 
in the fermionic case, by Frohlich et al."'^* and Jaksic 
et al.^""^^. The existence of a steady-state is rigorously 
proved by deriving explicit expressions for both the lesser 
Green's function and the current in the infinite time limit. 
Our method is exact in the interaction and needs neither 
Langreth rules, nor Dyson equations for Keldysh-Green's 
functions. It provides convergent expansions in terms of 
the interaction strength, i.e., it shows that the NESS is 
an analytic function of the interaction. 

Moreover, the proof of the steady-limit covers the two 
complementary but different transport scenarios: the 
partitioning approach of Caroli et al.^^ and the partition- 
free setting coined by Cini^^ . Let us briefly remind them 
here. In the partitioning approach the central region 
(sample) is coupled to biased leads at some initial in- 
stant. In contrast, the partition-free setup starts from a 
coupled and unbiased system, a bias being switched at 
t ^ to. Both settings have also received interest from the 
mathematical point of view^"^"^^'^^. 

The content of the paper is as follows: Section II sets 
the model and some notations needed for the partitioning 
and partition-free settings. We formulate our main re- 
sult, the steady state limit of the lesser Green's function, 
in Section III. The expression for the steady currents is 
derived in Section IV, along with the Landauer-Biittiker 
formula which holds if the Coulomb effects are considered 
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up to the 1st order term in the interaction strength. We 
outHne the proof of our results in Section V and conclude 
in Section VI with some general comments. 



II. 



MODEL AND NOTATIONS 



Our system consists of a finite sample S coupled to 
M semi- infinite leads labeled by 7. It is described by 
a discrete model: sites from the lead 7 are denoted by 
{*7}i>o and {m}mes arc the sites of the sample. Wc de- 
note by Hs and T-Ll the one particle Hilbcrt spaces of the 
sample and leads. The one particle Hilbert space of the 
compound system is % ~ Hl ® T~Ls and J- denotes the 
fermionic Fock space over %. The one particle Hamil- 
tonian of the noninteracting sample is hs, an arbitrary 
self-adjoint operator on T-Ls- In terms of the on-site cre- 
ation/annihilation operators on J-, the Hamiltonian and 
number operator of the sample are 



Hs 



771, n^S 



[m\hs\n) al^a^, = XI 



ni£S 



The semi-infinite leads are described by the Hamiltonian 
Hl = J2-y H-y and number operator Nl = Nj, where 



Hy=rLY. (al^a(,+i)^+ 4+i)."'J ' ^7 = H ^l". 



i>0 



i>0 



We also need the tunneling Hamiltonian Ht, some con- 
stant potential Hq applied to the leads and the electron- 
electron interaction in the sample V, 



M 



7—1 7—1 

V = - v{m,n)al^a^aia„. 



m,nGS 

Here, r is the coupling strength, G S the contact 
site of the sample with lead 7, v-y a constant potential, 
^ the interaction strength and v(rn,n) a pair potential. 
The differences ~ define the bias between the cor- 
responding leads. By convention all the vanish in the 
partitioning case, so that Hb — 0. For both partition- 
ing ('p')/partition-free ('pf') settings we assume that the 
switching (of the coupling to the leads Hr/of the bias 
Hb) happens suddenly at t ~ 0. A smooth switching 
can also be treated up to some technicalities^'^ , and does 
not influence the results. 

The Hamiltonians Hs, Hl, Ht, Hb as well as 



Hf) — Hs + Hl + Hi 



H — Hq + Hj 



act in the Fock space J- as the second quantized 
versions of single particle tight-binding Hamiltonians 
hs,hL,hT,hB,ho,h acting on H. Similar relations hold 
for the number operators Ns, N^, Nl and N = Ns + N^. 



We denote by Tl and Ts the subspaccs of T where 
Ns = and Nl = respectively. The full dynamics 
of the system is generated by the Hamiltonian 

K = H + V. 

We now introduce the initial state of the system as 
thermodynamic limit of states defined by density matri- 
ces on the sample coupled to finite leads of length A. 
Indicating this infrared cutoff by the superscript , we 
set 



(A) -I 



Tr (A){e^ 

for fl — [/ii, /iA/] and, for any leads observable Ol, 

A— ^00 

We recall that this state is characterized by the two-point 
function 

= ^a^{-i'f\Ff{h'f)\jt), 

and Wick's theorem. There, F^{e) = (1 + cf'i^~^|'-,))-^ 
denotes the Fermi-Dirac function of lead 7. 

In the 'p' setting, the initial state is the product state 
defined by 

{OlOs)p = {Ol)l.p (Os)s,p 

for any observables O^/s of the leads/sample. There, 

{Ol)l,p = {Ol)l,p, {Os)s..p = Tr:r,{ps,pOs}, 

where ps,p is an arbitrary density matrix on J-s- In the 
'pf ' case the leads and sample are coupled and have the 
same chemical potential pQ, i.e.. 



(0)pf = Um Tr^,A,{p(t^O}, 



A— >c>o 



where 



Ppf 



(A)n 



Tr^(A){e-''(^«+^i- +-^^^+^"^0^*"')} 
For later reference, we also define the leads state 

( ■ )L,pf = ( ■ ) L,p=[^io.....,tio]- 

The lesser Green's function is defined as 

where x,y are sites from either leads or sample, a^{t) = 
e'*^a^e~'*^ and □ stands for either 'p' or 'pf. This 
object plays a central role in the Keldysh approach and 
allows to compute both the particle density and the cur- 
rents. Our main concern being the existence of steady 
currents, wc are primarily interested in its large time be- 
haviour s — )■ 00 for constant t — s. 
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III. THE EXISTENCE OF THE STEADY-STATE 

For e v., we set 

aim = E('^i^) «^(i<^)) = E(^i<^) 

and let a#(|0)) denote either a{\(j))) or a^(|(/))). If g is 
a self-adjoint operator on Ti and Q denotes its second 
quantized version, then the well known identity 



a#(e'«|0)), 



holds. The interaction picture operators 

^f(t)=c-»MWe"^, 
satisfies the equation of motion 



(1) 



(2) 



a* +i 



-iuH 



Vc'''",A*{u)]du. (3) 



By Eq.(l) we have 

m.neS (4) 

xa^(e-'"''|m))a(e"'"'V))a'^(e"™''|n))a(e-™''|n)). 

Thus, the Dyson expansion of Af{t) obtained by iter- 
ation of Eq.(3) is a sum of iterated integrals involving 
monomials of the type 

A^(|0i),...,|0fc))=a#ni0i))---a*'=(l'/'fc)), (5) 

each being either \x) or e~'"''|m) for some m G S and 
u G [0,t] (so that, in particular < 1 for any 

unit vector j^f), |$) G T). Moreover, one easily sees that 
this expansion converges for any In fact, a careful study 
of this expansion shows that, under suitable assumptions, 
it remains convergent even for t = oo. 

Theorem III.l. Assume that the single particle Hamil- 
tonian h has neither eigenvalue nor real resonance. If the 
interaction strength ^ is small enough, then the limits 

lim A*{t), 

t— f oo 

exist. Moreover, a convergent expansion of Af in powers 
of the interaction strength ^ is obtained by setting t = oo 
in the Dyson expansion of A'^it). 

Remark. The first hypothesis of the previous theorem 
requires some comments. The spectrum of h^ + hs is 
continuous, filling the union of [v-y — 2tl, + 2tl]. If all 
the eigenvalues of the isolated sample hs are embedded 
in these bands and if the coupling to the leads r is weak 
enough then all these eigenvalues will generically turn 
into resonances of positive width. In such circumstances, 
the spectrum of h is continuous and coincides with that 



of h^ + hs. Moreover, one can show that for any a;, y in 
either the leads or the sample. 



\{x\e-''''\y)\dt < oo. 



(6) 



However, as the coupling t increases, some resonances 
may become real, cross a band boundary and turn into 
an eigenvalue of /i, invalidating (6). The first hypothesis 
in Theorem III.l is meant to ensure the validity of Eq.(6). 

Let w be an operator on such that K^luij-fA)! < 1 for 
all unit vectors |(/)), G "H. Replacing each term (5) in 
the Dyson expansion of Af hy M{w\(f>i), ... ,w\(j)k)) does 
not alter the convergence of this expansion. We denote by 
.A^[u'] the operator obtained from this modified Dyson 
expansion. 

We are now in position to state our main result (recall 
that /is = in the 'p' case): 

Theorem III. 2. Under the assumptions of Theorem 
III.l one has, for any s, 



(7) 



= i(At[e--(''-+"-)4]^Jc.;])i,n. 
There, w+ denotes the M0ller operator^ 

lim e'*^"'*''>L|<^), (8) 



where pl projects on the leads subspace % l • 

In Section V, we shall outline the proofs of Theorems 
III.l, III. 2. Complete mathematical details will be given 
clscwhere^'^. We conclude this section with several re- 
marks. 

1. Asymptotic completeness^ implies that 



41^-) 



lim e'*''°e-'*''|V') 



lim pLe'*"'^+''«)e-'*''|V'), 



(9) 



is unitary from % to TLl so that the object under the 
expectation on the RHS of Eq.(7) has a convergent ex- 
pansion as described above. Moreover, the expectation 
is w.r.t. the leads state ( • which does not depend on 
the interaction V , i.e., satisfies Wick's theorem. 
2. Eq.(7) implies right away that the expected particle 
number in the sample reaches a steady value in the long- 
time limit 



lim {Ns{t))a 

t^OO 



mes 



< 



(0). 



In fact, one can show that, under the assumptions of The- 
orem III.l, the system reaches a NESS ( • )□+ described 
by 



,#1 



■a*^)B+= }im{a*;{t)---al-{t))a 

= (.A#Mc.;].-.^#M4])i,H. 
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In the 'p' case, this NESS is independent on the initial 
state of the sample ps.p- Moreover, in the special case 
fl — [/io, ■ • • , Mo], ( ■ )p+ is the unique equilibrium state 
of the interacting system at inverse temperature /3 and 
chemical potential fxo- 

3. In the non-interacting case (V^ = 0), the NESS satisfies 
Wick's theorem with the two points function (aj^a^)0_|_ = 
{x\p+\y), where the one-particle density operator is 
given by^'^'^'* 



,t 



(10) 



IV. THE STEADY-STATE CURRENT 



The current operator of lead a 



where I lb assumes a Landauer-like form 
/o.LB=E r\F^{E)-F,{E))\r^^{E,OfdE (16) 

^ J-2TL 

with the transmittance^^ T^^{E,^) corresponding to a 
mean-field Hamiltonian /is,mf = + ^vmf where 

VMF = '"H^rn\m){m\ - ^ \m){n\ (17) 



mes 



and 



VH,m = ^ v(m, n){n\p+\n), vx.mn = v{m, n){n\p+\m), 



nes 



are Hartree and exchange terms, with the single-particle 
density operator p+ given by Eq.(lO). 



= -i[A-,7Vj = -i[iJ,7VJ 



is the second quantized version of the single-particle cur- 
rent ja ~ —i[h,pa], where pa projects on lead a. Its 
statistical average is^ 

(Ja(i))H =r(G<,„^oJt,0-G<„^„,„(t,i)). (11) 
We introduce the interaction picture current operator 



Jait) 



-itHAtK 



-itK itH 



(12) 



which is similar to Ax{t), Eq.(3) being replaced by 

Jcit) = Ja+i /*[e-™^ye'"^, Mu)]du. (13) 
Jo 

Using Eq.(7) with s = in Eq.(ll) we get 

= lim (J„(t))H = (14) 

where J'a = limt_j,oo Jait) is calculated by setting t = oo 
in the Dyson expansion generated by iteration of Eq.(13) 
and J7a[w^_] is obtained in the usual way from J^a- Com- 
paring the final formulas for the two cases 'p' and 'pf, 
one realizes that Ia,pi{v = 0) = la.pip = [mo,--,Mo])- 
Since = -i[H + V, N^] = -i[H^^^ + V, Ni^\ one has 

Ia,pf{v = 6) = - lim Tr:^,.,{plfi[H^^^+V,Ni^^]} = 0, 

due to the cyclicity of the trace and the fact that 
[p^J^\h'^^^ +V] ^ 0. Thus both currents vanish in the 
absence of bias. This fact cannot be seen from the inter- 
acting Meir-Wingreen formula^. 

The interaction efii'ects can be calculated pcrturbatively 
from Eq.(13). For the partitioning setting with identical 
leads having a hopping constant > one finds 



V. PROOFS 

We start by analysing the structure of Ax{t), 
foUowing^^. Iterating Eq.(3) and using Eq.(4), one ob- 
tains an infinite series involving iterated integrals of the 
multiple commutators 

[aim ("i) • ■ • fl"! (wi), [■ • ■ , [aln^iur] ■ ■ ■ an,(wr), a 



Xi i ? 



/a,p = /aXB+0(r)+0(Cr'^), 



(15) 



where «„(«) = a(e ™'^\ni)). By repeated use of the al- 
gebraic identity 

[bi • • ■6fe,ci • ■ -q] 
fe I 

= YY{-iy{h,Cj}ci ■ • -Cj-i^i • • -fi^ ■ ■ - bkCj+i ■■ - ci, 
1=1 j=i 

and of the canonical anti-commutation relations 
{at(|(/,)),a(|V'))} = im, am)} = we recast 

our expansion into the form 

Ax (t) = ax + 'YC dui • • • dur 

^y-^ JO<U,< - <Ul<t 

C^{ui,..,Ur;x)M'^{ui,..,Ur;x), (18) 

Ger^ 

where each is a finite set (of contraction diagrams). 
For each G G F^, is a monomial of type (5) and 

is a product of 'pairing factors' like (?/|e~™3^|?/') or 
(y|e-iK-«fc)'»|y')^ where y,y' eSU {x}. 

Our first assumption ensures that there exists a con- 
stant Cx such that 

max \{y\e-""'\y')\du<Cx. 

y,y'GSU{x} 

A delicate combinatorial analysis then shows that (see 
Theorem 1.1 in^^) 

^1^1'' / dui---du,. ^ \C^{ui,..,Ur;x)\ < oo, 

r>l ^0<Ur<-<iii<c<jQ^Y' 
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provided |C| < Aq = 2/ {27\S\^G^v), where \S\ is the num- 
ber of sites in the sample S and v = max„ \v{'n, m)\. 
Thus, the expansion (18) converges uniformly w.r.t. t G 
[0,oo]. In particular, setting i = cx) in Eq.(18) yields a 
convergent expansion of Ax (and taking adjoint gives an 
expansion for A\.). This proves Theorem III.l. 

To prove Theorem III. 2, we first notice that, according 
to Eq.(l), we get an expansion of 

by replacing each factor a^{\il))) of any monomial 

in Eq.(18) by a#(e-'*(''^+''«)e'(*-'')''|V')). Since Eq. (9) 

implies that 

lim a#(e-'*(''^+''«)e'(*-'')''|^)) = a#(e-'"(''^+''«'a;l|V')), 
one has 

hm e-"(^-+^-)af(t-s)e"(^-+^-) 

and hence, B{t) c-"(^^+^«)at (i - s)a^(t)e'*(^^+-f^«' 
satisfies 

lim B{i) = Al[c'''<-''^+''^^u;l]AJuji]. (19) 

Notice that since the range of is 7^^, the RHS of 
this identity is an observable of the leads. In the 'p' 
case, Hb = and the state ( • )p is invariant under the 
dynamics of Hl- It follows that 

lim (4(t-s)a,(t))p = hm (Z?(t))p 

which proves Theorem III. 2 in the 'p' case. 

To deal with the 'pf ' case, we invoke standard pertur- 
bation theory (see e.g.^°) to write 

There, ( • )d denotes the grand canonical ensemble for the 
decoupled dynamics Hq + V at inverse temperature /3 and 
chemical potential fiQ, i-C-, the product state 

and 

= / + rdri-.. 'drfci?T(Ti)..HT(Tfc), 

1.^1 "'0 JO 



where Ft(m) = e''(^"+^)i?Te-"('^"+^^ The state (•)d 
being invariant under the dynamics of + Hb , one has 

{Va^y{t- s)ax{t))d = (Ae(i))ci, 

where Vt = e-it{HL+HB)x)eiKHi^+HB) ^ follows from 
Eq.(19) that 

lini(P4(t-s)a,(t))d 

= lim (A4[e"'^(''-+'^-)c.|]^Jc.;])d. (21) 

Introducing the partial trace 

_ Tr^s{e-^(^s+y-MoWs)p| 
~ Trjrg{e-'3(^s+i^-MoJVs)} ' 

we observe that, for any observable of the leads, 

(AOL)d - (e-'*(^-+^-)2?ic"(«-+^«)0L)L,pf. 

Since h^+hB has continuous spectrum, the dynamics of 
Hl + Hb is mixing w.r.t. the state ( • ) L,pf (see e.g.-^^) so 
that 

lim {VtOL)d = fim (c-"(^-+^-)pie'*(^-+^-)OL)L.pf 

= {'DL)L,pf{OL)L,p{ = (^?)d(OL)L,pf. 

Applying this identity to the RHS of Eq.(21) and insert- 
ing the result into Eq.(20) proves Theorem III. 2 in the 
'pf case. 

VI. CONCLUSIONS 

We give sufficient conditions for the existence of a 
steady-state regime for open interacting systems, using 
an approach based on time-dependent scattering theory. 
While the approach is pcrturbative w.r.t. the interaction 
strength, the existence of the stationary regime is gen- 
eral in the sense that it does not rely on a particular 
approximation scheme for the Coulomb effects. To our 
best knowledge, the steady-state regime is not proved 
within the Keldsyh formalism: the steady-state value of 
the current is derived by assuming that the two-time GFs 
depend only on time difference^. We perform the ther- 
modynamic limit without making use of Langrcth rules 
and Dyson equations. In the Keldysh approach it is not 
clear how to perform the thermodynamic limit on the 
proper interaction self-energy. 

The smallness condition on ^ for the existence of the 
steady-state is also necessary. For example, if v{m, n) is 
diagonal then F is a single-particle operator (Hartree-like 
potential) which can create bound states for large values 
of leading thus to oscillations similar to the ones ob- 
served in^. However, we conjecture that if one performs 
the ergodic limit limT-s-oo ^ /q Jait)dt, a time-averaged 
steady-state could still be achieved even in the general 
interacting case. It would be interesting to investigate 
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the ergodic limit of the oscillating currents reported in 
Fig. 1 from^, which seem to support our conjecture. 

In the partitioning approach wc have shown that 
the steady-state quantities do not depend on the initial 
many-body configuration of the sample ps,p- Moreover, 
one can allow other switching procedures of the bias or 
of the lead-sample coupling (not just the sudden one), 
and the steady state remains unchanged (the complete 
proof will be given in Ref. 23). Let us mention that 
very recently^^ we have shown that when ^ is allowed to 
be arbitrarily large and the system is in the off-resonant 
regime in which h has eigenvalues situated very far from 
the continuous spectrum, the ergodic cotunneling current 
presents memory effects and depends on ps,p if £, ^ 0- In 
the non-interacting case, we still have independence on 

PS,p- 



Our results could be numerically implemented in both 
settings (partitioning and partition-free) and compared 
to the ones obtained from the Keldysh formalism. The 
second correction in Eq.(15) suggests significant differ- 
ences for strong coupling to the leads. 
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